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Abstract. We address the possible impact of New Physics on neutrino oscillation experiments. 
This can modify the neutrino production, propagation and/or detection, making the full 
cross section non-factorizable in general. Thus, for example, the neutrino flux may not be 
properly described assuming an unitary MNS matrix and/or neutrinos may propagate differently 
depending of their Dirac or Majorana character. Interestingly enough, present limits on New 
Physics still allow for observable effects at future neutrino experiments. 



1. Introduction 

In order to describe the possible impact of New Physics (NP) in the fuh neutrino oscillation 
process we must properly convolute the production, propagation and detection amplitudes 
because in contrast with the minimal Standard Model (SM), these subprocesses do not in 
general factorize for relativistic neutrinos. Comparing with the general expression for neutrino 
oscillations we can learn: (i) Under which conditions the production and detection processes can 
be described by pure states as defined by the Maki-Nakagawa-Sakata mixing matrix, or we must 
use the density matrix formalism, (ii) How good the factorization approximation is. (iii) If we can 
only consider NP effects in the neutrino propagation, neglecting departures from the SM in the 
production and detection processes, or viceversa. We find using a general effective Lagrangian 
including charged and neutral scalar and vector interactions for both neutrino chiralities, and 
constraining the new couplings to satisfy present bounds on muon decay and neutron beta decay, 
that future more precise oscillation experiments can be sensitive to NP. 



2. Main differences between the standard and NP approach 

Neutrinos with mass nii produced in the process la{^a) + A.{\a) — ^'i(A) -|- B{Xb) can be 
described by the density matrix (see Ref. [1] for notation): 

e%{p;\i-x\i') = j^^ J2 Anp,x-AB,XA,Xa)Ar{p,x'-AB,\AAa). (1) 

a 

This is equivalent to the usual description with pure Maki-Nakagawa-Sakata (MNS) states 
|i/q) = J2i=i U*^\l'i) if neutrinos have SM interactions and are relativistic, depending the 
difference between both descriptions on the NP strength. Present bounds on NP allow for 
effects at the few per cent level, thus at the reach of future, more precise neutrino oscillation 
experiments [l]. 



Then, the production density matrix evolves along the distance L = T to the detector with 
the Hamiltonian Tl describing the neutrino interaction with the surrounding medium, 

gfip, L = T = 0,)^ p%{p,L = T^O) = e"^^^ g^ip, L = T = 0) e*^^. (2) 

This generalizes the usual SM description [2|[3]. For instance, the different propagation through 
matter of Dirac and Majorana neutrinos in the presence of new right-handed interactions, which 
could be eventually observable, can be properly taken into account [4J. 

Finally, we can calculate the full cross section for the detection of the neutrino flavour /3 in 
the process + C ^ Ip + D from the density matrix at the detector gf)-. 

^ ' spins 

where E is the energy of the initial neutrino and the detection process amplitude 
[5]. This cross section factorizes for relativistic neutrinos within the SM, aa^p{E, L) = 
Pa-,f3{E, L) ap{E)^ where Pa^p{E^L) is the corresponding oscillation probability and crpiE) 
the detection cross section, but not in general. 

Hence, generically the number of neutrinos of flavour (3 produced per unit flux of neutrinos 
of flavour a and per scattering center is given by Eq. Obviously, the size of the new effects 
depends on the strength of the NP. If we try to be as much model independent as possible and 
only consider the experimental constraints derived from muon decay and neutron beta decay 
(see e.g. [6]), we find that future more precise neutrino oscillation experiments can be sensitive 
to NP. As an example we plot in the Figure {left) the possible departure from the SM prediction 
for the Vf. transition through the Earth for Dirac and Majorana neutrinos, and for L = 

732 km as a function of the incoming neutrino energy E [7] . Similarly for — > 1/^ in the Figure 
(right) but for propagation in vacuum, and for E = 10.5 GeV as a function of the baseline 
distance L. 
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